Suppression of Conductance in a Topological Insulator Nanostep Junction 
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We investigate quantum transport via surface states in a nanostep junction on the surface of a 3D 
topological insulator that involves two different side surfaces. We calculate the conductance across 
the junction within the scattering matrix formalism and find that as the bias voltage is increased, 
the conductance of the nanostep junction is suppressed by a universal factor of 1/3 compared to the 
conductance of a similar planar junction based on a single surface of a topological insulator. We 
also calculate and analyze the Fano factor of the nanostep junction and predict that the Fano factor 
saturates at 1/5, five times smaller than for a Poisson process. 

PACS numbers: 73.20.-r, 73.23.-b, 73.40.-c 
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Experimental demonstration of topological phases 
in both two-dimensional (such as HgTc) and three- 
dimensional (such as Bi2Se3) compounds with strong 
spin-orbit interaction [lHH has generated a plethora of 
interest in the physics community Q. These compounds 
are insulating in the bulk (since they have an energy 
gap between the conduction band and the valence band) 
but their surfaces support gapless topological excitations. 
These surface states are topologically protected against 
non- magnetic defects by time- reversal symmetry [5[ . 




FIG. 1: Schematics of the proposed junction, (a) The nano- 
step junction. As the name suggests the height L of the junc- 
tion is ~ 10 nm. (b) Similar junction, involving only one side 
surface of the topological insulator. See the text for further 
details. 

In the simplest case these low-energy excitations of a 
strong three-dimensional (3D) topological insulator can 
be described by a single Dirac cone at the center of the 
two-dimensional Brillouin zone (r point) [2, Si. Re- 
cently, naturally occurring defects such as step-like im- 
perfections 0,0] have been studied in the context of topo- 
logical insulators. Conductance through atomic scale 
step-defects has been measured demonstrating that the 
surface states are protected in spite of such an abrupt 



defect [fj . These atomic sized steps are usually modeled 
as delta function potential barriers Q . In this article we 
investigate quantum transport through a nanostep junc- 
tion involving two different side surfaces of a 3D topo- 
logical insulator. We predict that the conductance of 
the nanostep junction is suppressed in the large-energy 
limit by a universal factor of 1/3 as compared to the con- 
ductance of a similar junction based on a single surface 
of a 3D topological insulator (or a similar junction in 
graphene). Figure HJa) shows a schematic of the nano- 
step junction considered. The junction is divided into 
three regions. In region I (x < and 2 = 0) and in re- 
gion III (x > and z = L) the surface states lie in the 
x-y plane, while in region II {x — and < z < L) the 
surface states lie in the y-z plane. A dc bias voltage is ap- 
plied between region I and III and a top gate V s controls 
the carriers in region II. FigurcQJb) shows a schematic of 
the analogous planar junction, where in all three regions 
the surface states lie in the x-y plane and a top gate V p 
is applied to region II. The first junction [Fig. Ufa)] is 
referred to as a step junction, and the second junction 
[Fig. |TJb)] as a planar junction in the rest of this article. 
We use the conceptually transparent scattering matrix 
formalism to calculate the transport properties of these 
junctions For concreteness, we use typical parame- 
ter values of the 3D topological insulator Bi2Se3 when 
illustrating our results. 

The low-energy effective Hamiltonian for Bi2Se3 in the 
basis of four hybridized states of Se and Bi p z -orbitals 



denoted as 
written as 



,t), b2j,t>, \plt,i), Ip2 z ,|)} can be 



H(k) = e (k) 



M(k) 
Aik z 


A 2 k+ 



A x k z 
-M(k) 
A 2 k+ 






A 2 k_ 
M(k) 

—A\k z 



A 2 k- 


-Aik z 
-MQs) 



where k± = k x ± ik y , eo(k) = C + D\k\ + D 2 k+k_ 



(1) 



M(k) = M - B x k 2 z - B 2 k + k_, and fc + /j_ = k 2 x + k\ [l(J. 
Here f (I) stands for up (down) spin and + (— ) stands for 
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even (odd) parity. From this three-dimensional Hamilto- 
nian, there exists a straightforward procedure to obtain 
the effective Hamiltonian describing the surface states 
111 Il2| . The surface states in the x-y plane, for exam- 
ple, are obtained from the three-dimensional wavefunc- 
tions for these surface states (which are exponentially 
damped in the z-direction, with finite skin depth A) using 
Eq. ([I]), followed by imposing the boundary conditions of 
vanishing wavefunctions at the two boundaries (z = 
and z = L). For three-dimensional topological insulators 
(L ^> A) the surface states at the two boundaries are 
decoupled and the effective Hamiltonian describing the 
carriers in regions I and III is then given by (l2| : 



U xy = e xy + hv xy {a x k v - a v k x ), 
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where t xy = C + Q-M, hv\ 
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1 — -gi represents 



the Fermi velocity in the x-y plane and <r x , a y and a z 
denote the usual Pauli matrices. 

Analogously, we obtain the effective Hamiltonian de- 
scribing the carriers in region II as 
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-Aav/l - ^i- Solving 
the Hamiltonian we obtain the eigenstates in region 
I and region III as 
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FIG. 2: (Color online) Conductance as a function of energy 
for (a) the step junction [Eq. ([9|] and (b) the planar junc- 
tion [Eq. I|10p]. Different junction widths are plotted. In 
both plots, the solid (blue) line denotes L = 50 nm, the dot- 
dashed (green) line denotes L = 100 nm, and the dashed 
(red) line corresponds to L — 500 nm. Parameters used are 
Ai = 2.2 eVA, A 2 = 4.1 eVA, hv x F v = 4.065 eVA 



3.846 eVA, t% v = 0.03 eV, e y 



0.09 eV, eV s 



, hv v F z = 

-0.03 eV 



with corresponding energy eigenvalues given by e = €^ + 
hv x F \J k x + ky and tan(</>) = k y /k x . Similarly, in region 
II we obtain: 
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with corresponding energy eigenvalues e = £q Z + eV s + 

hvfJkl + {A\/Al)kl and tan( 7 ) = A 2 k y /(A 1 k z ). The 

+ (— ) labels of the wavefunction indicate right (left) trav- 
eling carriers in regions I and III, and downwards (up- 
wards) traveling carriers in region II. It should be noted 
that in general £g y 7^ £q Z , which implies that the Dirac 
cone describing the surface states in the y-z plane is 
shifted by an energy of £q = £q Z — £q with respect to 
the Dirac cone describing surface states in the x-y plane. 
Furthermore, the Dirac cone describing the excitations in 
region II has elliptic cross section (Ai 7^ A2). These fea- 
tures are in good agreement with recent electronic struc- 
ture calculations of similar systems [i~3| . 

Considering electrons incident from left to right, the 
total wavefunction in the different regions can be written 
as: 



vI/tt = a 1 ^ 
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in 



r*f if x < 0, z = 0, 
b^u if < z < L, x = 0, 
if x > 0, z = L. 



(6) 
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The reflection and transmission coefficients r and t can 
be obtained by imposing the boundary conditions un- 
der which the current normal to the boundary is con- 
served 14 - fl6| . We then hnd for the transmission proba- 
bility T = t*t of an electron incident on the step junction 
at a given angle of incidence (f>: 

cos 2 (</>) cos 2 (7) 
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tained by using conservation of energy and conservation 
of momentum along the y-direction. 

The zero-temperature conductance G s t op across the 
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FIG. 3: (Color online) Transmission (a) T step (0) [Eq. 0] for 
the step junction and (b) T p ianc(0) [Eq. (|10l) ] for the planar 
junction as a function of angle of incidence <f> for L — 100 nm 
and e = 0.25 eV. Parameters used are the same as in Fig. [2] 
The dashed red lines mark the envelopes of the transmission 
probabilities. 



step junction is then given by 
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stop 



Go 



it/2 



d(f)T stsp {4>) cos(0). 
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Here Go = ^p(c)kvFW, p(e) = e/(n(hvF) 2 ) denotes 
the density of states, W the sample width, and the in- 
tegration is over all angles of incidence <j>. Figure Eta) 
shows the conductance G s t ep [Eq. ©] function of 
energy e for different values of L. We observe that the 
conductance reaches a minimum at e = eo + eV s , arising 
from the difference in Dirac point energies in different 
planes. Beyond this point, the conductance first oscil- 
lates and then saturates at large energies. For k — > 1 
the conductance G s t cp reaches a limiting value of 2/3, 
independent of barrier width L. Below (see Eq. (fTTj) ) 
we discuss the behavior of the conductance close to this 
saturation point in more detail. 

The transmission probability of the analogous planar 
junction [see Fig. [Tfb)l with a top gate in the middle 
region is given by [17| : 



plane 
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cos 2 (0) cos 2 (7') 



(10) 

Here k' x represents the x-component of the momentum 



in region II, the energy dispersion is given by e 

,-1 



eV n 



-r hv x F v Jk'i + k 2 y and 7' = tan" 1 ^ ). Figure ^b) 
shows the conductance G p i anc as a function of energy e for 
different values of L. As before, the conductance reaches 
a minimum when e = eV p , and then increases to reach 
its saturation value G/Go = 1. For larger energies, the 
conductance of the step junction is thus suppressed by 
a factor of 1/3 compared with planar junctions. We at- 
tribute this suppression to the fact that the carriers in the 
step junction have to change their plane of propagation 
in region II [l8| . 

We now analyze in more detail the difference between 
the conductance of the step and the planar junction by 



comparing the denominators in Eqns. (J7J and (|T0l) . Fig- 
ure [3] shows the transmission probabilities T s t cp ((f>) and 
^ 1 piane(^ ) ) as a function of the angle of incidence 0. From 
Fig. (3Ja) we see that for the step junction there is a cut- 
off angle of incidence, which arises from the finite energy 
and velocity mismatch. This critical angle can be ex- 



pressed as 



'c,step 



sin 



We also see 



that the conductance of both junctions includes contri- 
butions from many resonant modes. Here, a resonant 
mode is defined as a mode with an angle of incidence for 
which the transmission T = 1. The various minima of 
these resonant modes form an envelope, as shown by the 
dashed (red) lines in Fig. [3] These envelope functions are 
obtained from the transmission expressions Eqns. (JT)) and 
(fTUfl by setting k z L = (2n+l)§ and k' x L = (2n+l)§, re- 
spectively, with n integer. When the number of resonant 
modes is large (ti>l)a lower bound for the integrated 
transmission is obtained by integrating over the envelope 
function. Under these conditions the conductance G s t ep 
of the step junction [Eq. (|9|] becomes 
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Similarly, we find for the planar junction 
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The suppression of the conductance 
Ggtcp in the vicinity of the saturation value thus depends 
on the shift of the Dirac point energies and the ratio of 
the Fermi velocities in region I and II [l9[. As final re- 
marks, we note that the effect of the elliptical dispersion 
in the middle region of the step junction can be incor- 
porated as an effectively wider barrier, L = and 
> 1, compared to the planar junction. The minimum 
barrier width needed to observe the suppression in the 
conductance described above is given by the condition 
for the existence of the first resonant mode in the junc- 
tion, L = ir/k z . In the case of Bi2Se3, this minimum 
width is L ~ 5 nm. For cleaved topological insulators, 
widths of L ~ 1 nm have been reported [20( . 

In the remaining part of this paper we investigate the 
Fano factor of the step junction, which is a measure for 
the noise suppression in the junction relative to Poisson 
noise Within the scattering matrix formalism the 
Fano factor F is defined as 



Jo /2 d^> cos(0) T(l-T) 



F = 



ll 12 d<p cos{<t>) T 



(13) 



where <j) is the angle of incidence of the carriers and T 
represents the transmission of the junction. Substituting 
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tion, namely the change of plane of propagation of the 
topological surface states. 

In conclusion, we have proposed and analyzed quan- 
tum transport through a nanostep junction on the surface 
of a 3D topological insulator in the ballistic limit. Our 
results show that the conductance in a nanostep junction 
is suppressed by up to a factor of 1/3 compared to similar 
junctions based on a single surface of a 3D topological in- 
sulator or graphcnc. Although the suppression depends 
on the ratio of the Fermi velocities and the difference in 
Dirac point energies of the different side surfaces of the 
3D topological insulator, the saturation values of the con- 
ductance and Fano factor themselves G s t cp — > 2/3 and 
F s tcp — > 1/5 are universal. We also predict oscillating 
behavior of the Fano factor around its saturation value. 



Experimental demonstration 22| of our predictions will 
provide further insight into the scattering mechanisms 
involved in topological insulator nanostep junctions. 
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FIG. 4: (Color online) The Fano factor [Eq. (JT3j>] for (a) 
the step junction and (b) the planar junction as a function 
of energy e. Different junction widths are plotted. In both 
plots, the solid (blue) line denotes L = 50 nm, the dot-dashed 
(green) line denotes L = 100 nm, and the dashed (red) line 
corresponds to L = 500 nm. Parameters used are the same 
as in Fig. [2] 



Eqns. ([7]) and (|10p into Eq. ([TBI and evaluating the in- 
tegrals in general leads to lengthy expressions. However, 
in the limit of a large number of resonant modes (n 3> 1) 
we find that the Fano factor of the step junction is given 

by [Hi 



R 



step 



12 
175' 



(14) 



where k as defined earlier. 

Figure @] shows the calculated F step and F p i anc for dif- 
ferent values of the junction width L. As the energy of 
the incident carriers increases, both Fano factors reach a 
maximum around e = 0.03 eV, where e = eo + eV s and 
e = eV p respectively. For higher energies the Fano factor 
of the step junction oscillates around its saturation value 
1/5, which is five times smaller than the Fano factor ex- 
pected for a Poisson process, F = 1. On the other hand, 
the Fano factor F p i ane for the planar junction vanishes as 
the energy of the incident carriers increases. This can be 
explained by noticing that the transmission T p i ane — > 1 
as e increases (see Fig. EJb)). This remarkable difference 
in the Fano factor for the two junctions again suggests 
that there exists an additional scattering mechanism for 
the step junction that does not exist in the planar junc- 
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